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Abstract 

We present in analytic form the matching conditions for the strong-coupling constant a^ s \^L) 



at the flavor thresholds to four loops in the modified minimal-subtraction scheme. Taking into 
account the present knowledge on the coefficient /?4 of the Callan-Symanzik beta function of quan- 
tum chromo-dynamics, we thus derive a five-loop formula for ai™-^^) together with appropriate 
£NJ ', relationships between the asymptotic scale parameters A&f' for different numbers of flavors rif. 
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The strong-coupling constant a s nf \fi) = g 2 /(4:ir), where g s is the gauge coupling of 
quantum chromo-dynamics (QCD), is a fundamental parameter of the standard model of 
elementary particle physics; its value ai (Mz) is listed among the constants of nature in the 
Review of Particle Physics Here, \i is the renormalization scale, and rif is the number 
of active quark flavors q, with mass m q <C /i. The /i dependence of ai n/ ^(/i) is controlled by 
the Callan-Symanzik beta function of QCD, 

^2 d «i W/) (/i) = (n/) / aj W/) (/x) 
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The calculation of the one- loop coefficient /? y about 33 years ago |2f has led to the discovery 
of asymptotic freedom and to the establishment of QCD as the theory of strong interactions, 
an achievement that was awarded by the 2004 Nobel Prize in Physics. In the class of schemes 
where the beta function is mass independent, which includes the minimal-subtraction (MS) 
schemes of dimensional regularization j^, (3^ f ^ and (3^^ |J| are universal. The results for 
& nf) and ^ n/) Q are available in the modified MS (MS) scheme Q- As for the 



term proportional to rij 
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where ( is Riemann's zeta function, was found in the large-n/ expansion |8j, while the residual 
terms, of O(n^) and below, are presently unknown. However, the latter were estimated by 
an educated guess, through weighted asymptotic Pade approximant predictions (WAPAP's), 
which are improved by including asymptotic corrections with respect to the usual Pade 
approximants and performing a weighted average over negative values of rif |9(. In the case 
of P^" f \ leaving aside the quartic Casimir terms, which appear there for the first time, the 
WAPAP's approximate the exact coefficients of ri"j with n = 0, 1, 2 amazingly well, at the 
one-percent level. One may thus expect that the WAPAP's for (3^^ work similarly well, 
except for the quartic Casimir terms, which cannot be predicted quite as reliably. For the 
reader's convenience, f3^^ (N = 0, . . . , 4) are listed for the rif values of practical interest in 
Table □ 

In MS-Hke ™* atio. S cne mes , the A P pe lqUlst -C_ decoupHng t Ke 0rem Q 
does not in general apply to quantities that do not represent physical observables, such as 



TABLE I: MS values of (3^^ for variable rif. (3^^ is estimated by WAPAP's with quartic Casimir 
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25 
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12 
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77 
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12 
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384 
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3456 

9769 
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65 
128 



445 s(o\ , 140599 
32 S W ' 4608 
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5184 \°J " r 373248 

11027 Afq\ _ 598391 
648 373248 

11237 a^o^ _ 63559 
576 4608 



162 
119 
107 
124 



beta functions or coupling constants, i.e., quarks with mass m q 3> /i do not automatically 
decouple. The standard procedure to circumvent this problem is to render decoupling explicit 
by using the language of effective field theory. As an idealized situation, consider QCD with 
rij = rif — 1 massless quark flavors and one heavy flavor h, with mass rrih ^> /i. Then, one 
constructs an effective nz-flavor theory by requiring consistency with the full n /-flavor theory 
at the heavy-quark threshold Jx = 0(rrih). This leads to a nontrivial matching condition 
between the couplings of the two theories. Although, a^irrih) = ai n/ ' ) (m/ l ) at leading and 
next-to-leading orders, this relationship does not generally hold at higher orders in the MS 
scheme, i.e., ai n/ ^(/i) starts to exhibit finite discontinuities at the flavor thresholds. If the /i 
evolution of ai n/ ^(/i) is to be performed at + 1 loops, i.e., with the highest coefficient in 
Eq. (^) being (3^ f \ then consistency requires that the matching conditions be implemented 
in terms of A^-loop formulae. Then, the residual [i dependence of physical observables will 
be of order A^ + 2. The QCD matching conditions at the flavor thresholds to two IllJ and 

n 

three |12j loops are known in analytical form; they are routinely used in the literature and 
even copied to the Review of Particle Physics jl(. Recently, the four- loop result was found, 
in semi-analytical form Q|. In fact, the most intricate four-loop tadpole master integrals 
involving one non- vanishing mass among the basic set that enters a ny s uch calculation could 
so far only be computed numerically, with limited precision 14, 2- ^ is the purpose of 



this Letter, to overcome this bottle-neck by presenting the four-loop matching condition for 
ai n/ ^(/i) entirely in terms of elementary transcendental numbers. This requires the analytic 
evaluation of the massive four-loop tadpole dia gram that is called X or T 91 in the recent 
literature [l|. Together with the results of Ref. (3], we thus enhance the knowledge of the 
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basic set of massive four-loop tadpole master integrals in analytic form. 

Prior to explaining the core of this analysis and presenting our analytic result for the 
four-loop matching condition for a« (//), we derive the five-loop formula for this coupling 
for fixed value of nf. In order to simplify the notation, we introduce the couplant a^ n f\fx) = 



ft 



(/j,)/tt and omit the labels /i and rif wherever confusion is impossible. Integrating 



Eq. (JU) leads to 
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where b^ = Pn/Po {N = 1,...,4), A is the so-called asymptotic scale parameter, and 
C is an arbitrary constant. The second equality in Eq. (J3J) is obtained by expanding the 
integrand. The conventional MS definition of A, which we adopt, corresponds to choosing 

. Iteratively solving Eq. Q yields, with L = ln(/i 2 /A 2 ), 
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The particular choice of C [?1 Il2| in Eq. (j3J) is predicated on the grounds that it suppresses 
the appearance of a term proportional to (const. /L 2 ) in Eq. (@J). 

We now turn to the analytic evaluation of the four-loop matching condition for qg (//) at 
the flavor thresholds. The underlying formalism was comprehensively explained in Refs. 11, 
121 ]. and most of the technical issues related to its application at four loops were already 
discussed in Ref. For lack of space, we thus concentrate here on the missing link of 

this analysis beyond the scope of Ref. [13j], namely the analytic evaluation of the massive 
four- loop tadpole diagram X Q , which is depicted in Fig. Ufa). This task may be simplified 
by noticing that X does not represent a master integral, but may be reduced to simpler 





FIG. 1: Four- loop tadpole diagrams (a) X$ and (b) Jo- Dashed and solid lines represent massless 
and massive propagators; a dot on a line duplicates that propagator. 

integrals with less lines, all of which are analytically known some for a short time only 
jlfil |. except for the one (Jo) shown in Fig.^b). The integral Jq is finite, and the coefficients 
of its expansion in e, where D = 4 — 2e is the dimensionality of space-time, have only one 
level of trancendentality [3], i.e., they contain poly-logarithms Li^ and zeta functions ((k) 
with the same value of k. These properties reduce the number of terms and thus simplify 
the calculation. In order to evaluate Jq, we temporarily introduce an artificial mass splitting 
among the four massive lines in Fig. ^b), by assigning the mass m to any two of them 
and the mass M to the other two. We then perform an expansion in the ratio x = m 2 /M 2 
using the large-mass expansion technique and recover the complete series in x as explained 
in Ref. Q|. Through 0(e 2 ), we have 

CO f 

(1 - e)(m 2 M 2 ) 2e m 2 J = ^x n I - 

n=l ^ 

C(3)-5 3 , 
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where we introduced the short-hand notation S a = 2 a 2 S a {2n — 1) — S a (n — 1), with S a (n) = 
Y^j=\j~ a being harmonic sums, and omitted irrelevant terms involving lnx. We then put 
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x = 1 and exploit the identities 
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x [5 r 6 (Z)+2S r _ fc (Z)], 
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at trancendentality levels k = a + c and = a + 6 + c, respectively. The sums with k < 5 
may be found in Ref. while those with k = 5 may be obtained from there through 
integration, 

f(l)x l _ f x dy^ f(l)y l 



^ i c+l Jo y ^ 
i=i Jl) y i=i 



where f(l) = S± a (l), . . .. After some algebra, we find an analytic expression for J and hence 
also for X , 
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X = -318C(4) In 2 + — C(5) - 486 5 + 0(e). 



Here and in the following, we use the constants 

( , x In 2 2 
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If we measure the matching scale \i in units of the MS mass (/i), our result for the 
ratio of a! = aS n ^ (p) to a = a^ nf ' (fx) reads 
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FIG. 2: jx dependence of cts(Mz) from iV-loop evolution and (N — l)-loop matching, with N = 1 
(dotted), 2 (short-dashed), 3 (dot-dashed), 4 (long-dashed), and 5 (solid). 



where i = In [li 2 /m\ (fj,)] and 
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The counterpart of Eq. fTU|) in the on-shell scheme of mass renormalization may be obtained 
by substituting the three- loop relationship between rrih(Li) and the pole mass Mh jlfll ]. 

Going to higher orders, one expects, on general grounds, that the relationship between 
ai"'^//) and a™ (fx), where // <C Ji <C li, becomes insensitive to the choice of /2 as long 
as Ji = 0{my l ). This has been checked in Ref. 12| for four-loop evolution in connection 
with three-loop matching. Armed with our new results, we are in a position to explore 
the situation at the next order. As an example, we consider the crossing of the bottom- 
quark threshold. In particular, we wish to study how the Jjl dependence of the relationship 
between a^\M T ) and ai 5 \iMz) is reduced as we implement five-loop evolution with four- 
loop matching. Our procedure is as follows. We first calculate af 1 (ji) with Eq. (JU) by 
imposing the condition ai 4) (M T ) = 0.34 , then obtain ai 5 \/x) from the on-shell version 



7 



of Eq. (dm with M b = 4.85 GeV [if], and finally compute af\M z ) with Eq. ®. For 
consistency, AT-loop evolution must be accompanied by (N — l)-loop matching, i.e., if we 
omit terms of 0(1/ L N+l ) in Eq. (j3j), we need to discard those of 0(a ) in Eq. (jlOJ) at the 
same time. In Fig. El the variation of ai 5 \Mz) with p./ Mb is displayed for the various levels 
of accuracy, ranging from one-loop to five-loop evolution. For illustration, \x is varied rather 
extremely, by almost two orders of magnitude. While the leading-order result exhibits a 
strong logarithmic behavior, the analysis is gradually getting more stable as we go to higher 
orders. The five-loop curve is almost flat. Besides the fi dependence of af' (M z ), also its 
absolute normalization is significantly affected by the higher orders. At the central scale 
ft = Mf,, we encounter an alternating convergence behavior. 

As we have learned from Fig. in higher orders, the actual value of fl does not matter 
as long as it is comparable to the heavy-quark mass. In the context of Eq. ()1U|). the choice 
/2 = where Hh = m h (fi h ) is the renormalization-group (RG) invariant MS mass, is 
particularly convenient, since it eliminates the RG logarithm £. With this convention, we 
obtain from Eqs. (jSJ), ®, and (jlOJ) a simple relationship between A' = A*™') and A = A( n t\ 
viz 

# In ^ = (# - A,) I + (b[ - h) In I - b[ In ^ (12) 

+ in h ( b 'i ~ & i) ln/ + b 'i ~ b i- b 2 + b 2 + c 2 ] 
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x^|- C3 )-^-i(^ fo4 ) + c 4 } + o(^ 
where I = ln(/i 2 /A 2 ). The C(l// 3 ) term of Eq. (|T2^) is new. Equation (fT2^) represents a 
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closed four-loop formula for in terms of A( n ^ and /jl^. For consistency, it should be 
used in connection with the five-loop expression (J3J) for a^ f \^) with the understanding that 
the underlying flavor thresholds are fixed at p, = /Jh- The inverse relation that gives A^' 
as a function of A*™') and fXh emerges from Eq. (|T2j) via the substitutions A <-> A'; /3 N <-> (3' N 
for JV = 0, . . . , 4; and c N -> -cyv for V = 2, 3, 4. 

In conclusion, we have extended the standard description of the strong-coupling constant 
in the MS renormalization scheme to include five-loop evolution and four-loop matching 
at the flavor thresholds. These results will be indispensable in order to relate the QCD 
predictions for different observables at next-to-next-to-next-to-next-to-leading order. 
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